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Abstract 



A measurement of the branching ratio for the rare decay mode Bs A*^/^^ at 
the Tevatron is an opportunity to test various supersymmetric scenarios. We 
investigate the prospects for studying this mode in Run II and estimate that 
CDF would be sensitive to this decay for a branching ratio > 1.2 x 10^*^ with 
15 fb~^(or, if a similar analysis holds for DO, > 6.5 x 10^^ for the combined 
data). We calculate the branching ratio in minimal supergravity (mSUGRA) 
parameter space, and find that tan /? > 30 can be probed. (This mSUGRA 
parameter space cannot be probed by direct production of SUSY particles at 
Run II.) Including other experimental constraints on the mSUGRA param- 
eter space, one finds that CDF Bs fJ-^fJ-' measurements would be able to 
cover the full mSUGRA parameter space for tan /3 = 50 if the muon — 2 
anomaly exceeds ~ 11 x lO"^'', and about half the allowed parameter space 
for tan/3 = 40. A large branching ratio > 7(14) x 10~^ (feasible with only 2 
fb~^) would be sufficient to exclude the mSUGRA model for tan/3 < 50(55). 
Dark matter neutralino-proton detection cross sections are examined in the 
allowed region, and should be large enough to be accessible to future planned 
experiments. Combined measurements of Bs — > the Higgs mass nih 

and the muon — 2 anomaly would be sufficient to determine the fi > 
mSUGRA parameters (or show the model is inconsistent with the data). We 
also briefly discuss the Bs decay in R parity violating models. There, 

for some models, the branching ratio can be large enough to be detected even 
for small tan (3 and large Tni/2- 
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1. INTRODUCTION. Supersymmetric (SUSY) extensions of the Standard Model (SM) 
represent a natural candidate for the new physics expected to occur in the TeV energy 
domain. One of the difficulties in determining predictions of such models lies in the large 
number of new parameters the theory implies. Thus the most general low energy model, the 
Minimal Supersymmetric Standard Model (MSSM), has over 100 free parameters, and even 
the minimal model based on supergravity grand unification at the scale Mq = 2 x 10^^ GeV 
(mSUGRA) possesses four new parameters and one sign in addition to the SM parameters. 
Fortunately SUSY models apply to a large number of different accelerator and cosmological 
phenomena, and a great deal of effort has been involved in recent years to to use this data 
to limit the parameter space. Part of the difficulty in doing this resides in the success of the 
model in not disturbing the excellent agreement of the precision tests of the SM [|l|] due to the 
SUSY decoupling theorems which suppress SUSY contributions at low energies. Historically, 
the absence of flavor changing neutral currents at the tree level played an important role 
in the construction of the SM. They represent therefore an important class of phenomena 
that might show the presence of new physics, since the SM and the SUSY contributions 
contribute flrst at the loop level with comparable size. Thus the decay b sj has been a 
powerful tool in limiting the SUSY parameter space. In this paper, we consider the decay 
Bg — > fi~^fi~ within the framework of mSUGRA models and R parity violating models. This 
process is particularly interesting for several reasons: The SM branching ratio is quite small, 

1. e., Br[Bs yU^/^^]sM = 3.5 x 10^^ [||. The SUSY contribution [^|-|^ has terms that grow 
as tan^/5 (where tan/5 =< H2 > / < Hi > and Hi 2 are the two SUSY Higgs bosons) and 
thus can become quite large for large tan (3. Finally, as we shall show below, the two collider 
detectors at the Tevatron will be sensitive to this decay for tan/3> 30 in Run II. 

2. MSUGRA MODELS. We consider flrst mSUGRA models with R parity invariance, 
and combine the existing accelerator constraints on the parameter space {e.g., the light Higgs 
mass rrih bound, b s'j etc.), the cosmological dark matter constraints, the possible muon 
magnetic moment anomaly, with what might be expected from the Tevatron measurement of 
the Bs —>■ fJ'~^fJ'~ decay. The combined constraints can signiflcantly limit the SUSY parameter 
space, and thus allow better predictions as to what the models predict at the LHC. 

The mSUGRA model PJTO| depends on four parameters which we take to be the fol- 
lowing: rriQ (the universal scalar mass at Mq), mi/2 (the universal gaugino mass at Mq), 
Aq (the universal cubic soft breaking mass at Mq), and tan/? (at the electroweak scale). 
In addition, the sign of /i (the Higgs mixing parameter appearing in the superpoten- 
tial as fiHiH2) is arbitrary. We consider the parameter range of mQ,mi/2 < 1 TeV 
3 < tan/3 < 55, and \Aq\ < 4mi/2- We take a 2cr bound on the b ^ s'y decay |Tl 



of 1.8 X 10"^ < Br[b — > 57] < 4.5 x 10""^, and the LEP bound on the light Higgs of 
nih > 113.5 GeV ||12[. Since there is still a (2-3) GeV uncertainty in the theoretical calcu- 
lation of TUh, we will (conservatively) interpret this to mean m/^ (theory) > 111 GeV. For 
tan/3 > 45, results are sensitive to the precise values of nih and rrit. We assume here that 
nT-blmb) = 4.25 GeV and rritipole) = 175 GeV, and include the large tan/? corrections to 
the 6-quark Yukawa coupling constant |]TB|]. We assume that any muon g — 2 deviation from 
the SM is due to SUSY [|l4| , pr5| . The experimental deviation has now been reduced to a 



1.6 a effect [|T6],|T7|, and we take here this deviation to be greater than la below the current 



central value i.e., a^^^^ > 11 x 10 We will show however, what would happen if this 
deviation were to change with the new BNL E821 data currently being analyzed. For mod- 
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els with R parity invariance, the lightest neutralino, Xi is the dark matter particle, and we 
require that 0.07 < Q-oh"^ < 0.21, in accord with current CMB and other astronomical data 

Al 



|T^. All stau- neutralino co- annihilation channels ]TP|-pT[] are included in the relic density 
calculations. 

3. Bg — > DECAY. The branching ratio for Bg fi^fi~ is given in |Q which we 

write in the form 
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where /^^ is the Eg decay constant, Mb is the B meson mass, tb is the mean life and mi is 
the mass of lepton. Cs, C5, Cp, C'p include the SUSY loop contributions due to diagrams 
involving the particles such as stop, chargino, sneutrino, Higgs etc.. For large tan/?, the 
amplitude has terms that grow like tan^ (3 as can be seen in the example of Figure |T]. Thus 
at large tan/3, the dominant contribution to Cs is given approximately by 
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where mr are the two stop masses, and 67 is the rotation angle to diagonalize the stop 

tl,2 t 

mass matrix. We need to muliply the above expression by 1/(1 + ebY to include the SUSY 



QCD corrections, e?, is proportional to /i tan/3 p2[. We have Cp = —Cs,Cg = {■ms/mb)C, 



and C'p = —{ms/mb)Cp. The operators are given in ref. [Q. In our numerical calculation, 



we use all complete one loop contributions to the branching ratio. 

4. DETECTION OF at TEVATRON. We consider now the possibility 

of detecting the decay Bs —* fi^fi^ by the CDF and DO detectors at the Tevatron in Run II. 
Both detectors have been upgraded with excellent tracking and muon detector systems 0. 
The dimuon trigger is the key to collect the Bg fi^fi~ decays. 

In order to estimate the limits on Br[Bs — > /i^/U^] detection, we use the 95% CL. limit 
on Br[Bs — > fi~^fi~] published by CDF Thus our discussion is based on the CDF 

detector, although both CDF and DO detectors should have a similar perfromance. 

In the Run I analysis, CDF observed one candidate that was consistent with Bg fJ'^f^~ 
with an estimate of 0.9 backound (BG) events in 98 pb^^ [^. The primary Run-I selection 
variables and cut values were that at least one muon track with cr = L^.yMp/PjJ^ > 100 fim, 
I = P^^'/[P^'' + EPt] > 0.75 for the muon pair, and A$ < 0.1 rad. Here, L^y is the 
transverse decay length; pjt^ is the transverse momentum of the dimuon system. UPp is the 
scalar sum of the transverse momenta of all tracks, excluding the muon candidates, within 
a cone of AR e 



[ArjY + = 1 around the monetum vector of the muon pair. The 

z coordinate of each track along the beam line |^ must be within 5 cm of the primary 
vertex. A$ is an opening azimuthal angle between Pp^ and the vector pointing from the 
primary vertex to the secondary vertex (the reconstructed 5-meson decay position). As a 
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conservative estimate, CDF took the one event as signal to calculate 95% C.L. limit of signal 
events (Nf^"^" =5.06 events |^) and had set a limit of Br < 2.6 x 10~^. In the analysis, the 
selection efficiency (e) for signal events and the rejection power (7^) for background events 
(pass a baseline selection [^) are estimated to be ei = 0.45 and TZi = 440 by using a sample 
of like-sign dimuon events (5 < M^^ < 6 GeV/c^). 

A dimuon trigger in Ref. will improve the acceptance for signal events by a factor of 
2.8. The trigger will soon be tested using the Run Ila (2 fb~^) data. This will allow us to 
modify the trigger design for the higher luminosity expected in Run lib (15 fb~^). In our 
analysis, we assume that the dimuon trigger can be designed by maintaing the acceptance 
for signal events. We expect to improve the acceptance for signal events by a factor of 2.8 
0. If we assume the factor 2.8 to be the same for BG events, then we would observe 51 
(386) events in 2 (15) fb"^ with the same cuts as in Run I. Therefore, CDF has to require 
a set of tighter cuts to obtain the best possible upper limit. 

We first need to understand the background contents and expected improvement by the 
new Run II detector. Two types of backgrounds must be taken into account: (i) non-6 
backgrounds comming from the primary vertex; (ii) b background events, such as the gluon- 
spliting bb events. 

The most important feature for B decays is the displaced vertex. One way to reduce 
prompt background is to require a minimum decay length L^y. However, two tracks can 
appear to form a secondary vertex if one of two tracks originates from the primary vertex 
and the other has an impact parameter (d). Therefore, the requirement of a minimal impact 
parameter of individual tracks can further clean up the sample. It has been shown for 
example, in the Run-I analysis for B^ — >■ K^*ix^ ii~ events |]2^, that a tight impact parameter 
cut on significance for individual track {d/ad > 2) significantly improve the background 
rejection even with L^y > 100 fim. One has then e ~ ei x 0.43 and TZ ~ TZi x 190. 

We indeed observed the similar numbers of OS and LS events in the Run-I final selection 
for Bg —>■ fi~^fi~ events, while we should expect more OS events than LS events if the 
cuts were tight enough for non-fe background events. This indicates the final dimuon event 
candidates in Run I analysis are actually dominated by the non-6 background sources. Thus 
a higher track impact parameter is neccessary to reduce the non-6 backgrounds. We would 
expect larger reduction with good efficiency even after the L^i-y cut. In other words, we 
would have been able to set better limits if we had included the impact parameter for the 
cut optimization. The silicon vertex detector (SVX-II) will provide us much better reduction 
for the non-6 background than Run I. The non-6 background will not be a problem in Run 
II. 

In Run II, the most severe background will be the two muons from gluon-spliting 66 
events. Since both particles are 6 quarks, the impact parameter does not help. Both 6 and 
6 also go in the same direction, so that cut on L^y does not help either. However, A$ is 
still usefuU to remove the background events. Furthermore, in Run II, we can use AO in 
r-z view since we have z-strips in SVX-II. 

There is some room to improve the isolation cut. We can form a new isolation parameter 
by only using the tracks with large impact parameter. This new isolation cut will work to 
reject the 66 rather than non-6 background. Furthermore, we can search for tracks with 
large impact parameter on the opposite side of the dimuon candidates to make sure that the 
6 and 6 go to the opposite side. 
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Therefore, CDF could improve the BG rejection by a factor of 200-400 with further 
reduction of the signal efficiency by a factor of 2-3. Based on these facts, we now consider 
two cases to evaluate Run II limits as a function of luminosity. 

In the first case (Case A), we naively assume new tighter cuts in Run II, described 
above, will gain additional BG rejection power of 450 for additonal efficiency of 0.45, or 
7^2 = 450°-^^/^2_ This gives us 



£2 1 



ei l + log(7^2/7^l)/log(450) 



(3) 



If we could optimize the BG rejection in Run Ila (2 fb^^) to be Tl2 ~ blTZi with £2 ~ 0.61ei 
(from Eq. then we would expect one BG event in 2 fb^^. Thus, with an assumption of 
the same size of the total systematic uncertainty in Run II as in Run I, we can extrapolate 
the 95% C.L. limit to be Br < 7.7 x 10^* for 2 fb'^ using A^f 

In the second case (Case B), we simply assume the Run-II background rejection could be 
improved (without loosing the signal efficiency) to keep the expected BG events in 2 fb~^ at 
the level of Run I (i.e., 0.9 events). If we would observe one event in 2 fb~^, then we could 
set the limits by scaling the Run-I Br limit down by the luminosity (2000 pb^^/98 pb^^) 
and the acceptance by (2.8/1.0). Thus we obtain Br < 4.6 x 10^^. This would certainly 
be the optimistic scenario, but it would be a goal of this analysis in Run Ila. Here, the 
systematic uncertainty in Run II is assummed to be the same as in Run I. 

We repeat the same arguement for different luminosities. Figure ^ shows 95% C.L. limits 
on Br[Bs /U^/W^] at CDF in Run II as a function of integrated luminosity for Cases A 
and B. For 15 fb^^ in case A, CDF would be sensitive to Br > 1.2 x 10~^ and the combined 
CDF and DO data (30 fb^^) would be sensitive to Br > 6.5 x 10"^. 

5. RESULTS FOR MSUGRA. We examine first the parameter region for the mSUGRA 
model that would be accessible to CDF at Run II with 15 fb~^ of data. Figure ^ shows 
the Br[Bs as a function of mi/2 for Aq = 0, mo= 300 GeV. One sees that with 

a sensitivity of Br[Bs > 1.2 x 10^^ for 15 fb~ , the Tevatron Run II can probe 

the Bg decay for tan /3 > 30. Further, a search for this decay would sample much 

higher regions of mi/2 than a direct search at Run II for SUSY particles which is restricted 



to mi/2 < 250 GeV |2g]. As mo increases, the branching ratio goes down. However, this 
dependence becomes less significant for large mi/2, where mo as large as 800 GeV can be 
sampled for large mi/2- 

In Figure H the contours of Br[Bs — > are plotted in the mo-mi/2 plane for tan/? 

= 50, Aq = 0. We combine now this result with the other experimental constraints. Thus 
the shaded region to the left is ruled out by the 5 — >• S7 constraint, and the shaded region 
on the right hand side is disallowed if a^^*^^^ > 11 x 10^^°. The narrow shaded band in the 
middle is allowed by the dark matter constraint. We note that independent of whether the 
astronomically observed dark matter is SUSY in origin, the dark matter allowed region for 
mSUGRA cannot significantly deviate from this shaded region, for below the narrow shaded 
band, the stau would be lighter than the neutralino (leading to charged dark matter), while 
above the band, mSUGRA would predict more neutralino dark matter than is observed. 

Using our estimate that Br > 1.2 x 10^^ can be observed with 15 fb~^, we see that 
almost the entire parameter space allowed by the constraint can be probed in Run II 
for tan/5 = 50. Note that an observed Br[Bs t^~^f^~] > 7 x 10~^, possible with only 2 
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fb~^(see Figure |), would be sufficient to rule out the mSUGRA model for tan/? < 50. In 
Figure ^ we also show the expected dark matter detector cross section for Milky Way dark 
matter (the short solid lines). They are of a size that can be observed by future planned 
dark matter detectors such as GENIUS, Cryoarray, ZEPLIN IV and CUORE. In Figure |^ 
we plot the same information for tan (5 = 40, Aq = 0. We see here about half the parameter 
space can be scanned by the CDF detector (and the whole parameter space if a similar 
analysis holds for the DO detector). See Figure ^. We note, further, that if Aq = 0, a. 
simultaneous measurement of both Br[Bs — > fJ'~^f^~] and would essentially determine the 
mSUGRA parameters, as the rriQ allowed region at fixed mi/2 is very narrow due to the dark 
matter constraint. The effect of varying Aq is shown in Figure where the allowed region 
for Aq = — 2mi/2, tan/5 = 40 in the mQ-mi/2 plane is plotted. The effect is to tilt (and 
narrow) the allowed dark matter band. The entire allowed parameter space can again be 
probed. 

In Figures |^, ^ and |^, we have also drawn lines for various light Higgs masses (vertical 
dotted lines). A measurement of Bs fJ'~^f^~ , cl^j, and m/j would then effectively determine 
the parameters of mSUGRA for /i > by requiring that they intersect with the dark matter 
allowed band at a point. (If no choice of parameters allowed this, mSUGRA would be ruled 
out.) The Tevatron Run II should be able to either rule out a Higgs mass or give evidence 
for its existence at the 3 a level over the entire allowed mass range of SUSY light Higgs 
masses. Alternatively, the LHC's determination of rrih or the gluino mass (to determine 
mi/2) would fix the parameters of mSUGRA. 

The BNL E821 experiment should shortly have a more accurate value for a^, with errors 
reduced by a factor of two or more. We note from the above figures the importance this 
result might have. Thus if increases, the bound moves downward, encroaching further 
on the allowed part of the parameter space, and a value of > 50 x 10^^*^ would eliminate 
the mSUGRA model ||2^. However, if decreases significantly (but is still positive), the 



mSUGRA model would predict a heavy SUSY particle spectra closer to the TeV region, 
having significant effects on accelerator and dark matter detection physics. An accurate 
determination of corresponds to a line from upper left to lower right (or more precisely 
a band when errors are included) running parallel to the < 11 x 10^^*^ boundary, and 
cutting through the allowed dark matter band which runs from lower left to upper right. 
Thus these two experiments are complementary for determining the mSUGRA parameters. 
6. R PARITY VIOLATING MODELS. We consider briefly here the case when R 
parity is broken. The general renormalizable R parity violating superpotential has the form 

= K,L,H2 + \,uUL,El + K^,UQ,Dl + K.iJJtD^l (4) 

The first three terms are lepton number violating, and the last term is baryon number 
violating. For the Bs — >■ /i"*"/!" decay, we need to consider lepton violating terms [^, and so 



we set A" to zero (to prevent rapid proton decay). Unlike the R parity conserving models, the 
SUSY contribution to Bg — > can now occur at the tree level, and so can be considerably 

larger. For example, if the A and A' terms are present in the theory, we can have 

ri ^j23^i22 ^j23^i22 (r\ 

In Figure |^, we plot the Br as a function of mi/2 in the R parity violating scenarios. We 
have chosen A^23 — ^i22 = 0.02 and tan/3 = 10. These small couplings are not restricted by 
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any another physical process |2^. We see that in this case smaller values of tan/3 can be 



probed and the Br is large compared to the R parity conserving case. The contributions of 
Cs,p to the newly discovered decay mode Bg Kfifi are small. For example, for mi/2 = 
300 GeV and mo = 300 GeV, 5r[S, K/i/i] gets a contribution of 3.30 x 10"^. (This 
calculation, however, involves a large QCD uncertainty.) The observed branching fraction 
for this mode is (0.99t[J:^l[!:J^) x lO^^ g. 

In contrast to the R conserving case, the R parity violating scenarios are not infested 
with tau's and therefore has the potential to be observed at Run II directly. For example, let 
us consider the production of xtx2 assume that Xi22, A^23 present. In this case 

we can have several interesting final states such as 61 + 2jets+missing Et (from neutrino) 
and 4/ + 4jet+missing Et from the production of xf-X2- 

7. CONCLUSIONS. We have investigated the prospects for studying the rare decay 
mode Bg fi^^^ by the CDF detector in Run II at the Tevatron. We have analyzed the 
background for this process and find that a Br > 1.2 x 10~^ for 15 fb~^ can be probed at Run 
II. (This is nearly a factor of 100 improvement over the Run I bound.) In this background 
analysis, we have not only extrapolated the Run I limit, but also introduced new selection 
cuts to reduce the background. The improvement of the Run II detector has also been taken 
into account. If a similar analysis can be performed for the DO detector, the combined data 
would increase the sensitivity for detecting this decay mode to Br > 6.5 x 10~^. The LHC 
should be able to reach the sensitivity at the level of SM branching ratio |^ . 

The Bs fi^fi^ is an important process for SUSY searches for new physics, as the Stan- 
dard Model prediction of the branching ratio is quite small (3.5 x 10~^), and the SUSY 
contribution increases for large tan (3 as tan^ f3. For the mSUGRA model, the above sensi- 
tivity implies that Run II could probe a region of parameter space for tan /? > 30, a region 
which could not be probed by a direct search at Run II. We have combined the expecta- 
tions for Bs fJ'^fJ'^ for mSUGRA with other experimental bounds on the parameter space. 
Thus a large branching ratio, i.e., > 7 x 10~^(14 x 10~^) would be sufficient to eliminate 
the mSUGRA model for tan/5 < 50(55). A measurement of Bg — > /^"""/i" , the muon g — 2 
anomaly and the light Higgs mass combined with the astronomical bounds on cold dark 
matter would essentially determine the mSUGRA model, allowing predictions of all other 
sparticle masses and cold dark matter neutralino-proton cross sections. If the muon g — 2 
anomaly remains positive these cross sections are ^10^^ pb, and should then be accessible 
to future planned detectors such as GENIUS, Cryoarray, ZEPLIN IV and CUORE. 

A simple model of R parity violation was also considered, and here the Br[Bs Ai+yU^] 
can be large for both small and large tan/?. 
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s 

FIG. 1. Example of diagram contributing to Bs — > fJ-'^fJ-' with leading contribution of tan^ f3. 
The H and A are the heavy CP even and CP odd neutral Higgs bosons, = Ij 2) are charginos 

and ti{i = 1, 2) are the stop bosons. Heavy marked vertices each contain a factor of tan p. 
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Integrated luminosity (fb~') 
FIG. 2. Illustrated 95% C.L. limits on the branching ratio for Bs — > at CDF in Run II 

as a function of integrated luminosity. Solid (Case A) and dashed (Case B) curves are based on 
different assumptions on the signal selection efficiency and the background rejection power. See 
the text for details. 
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FIG. 3. Branching ratio for Bg — > Ai"*"/^" as a function of for various tan /J values in 

mSUGRA models. Other mSUGRA parameters are fixed to be tjiq = 300 GeV, = and fj, > 0. 
Dashed and dash-dotted lines are to indicate the models are excluded via constraints on Br [b 57] 
and nif > m~^o, respectively. 
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FIG. 4. Branching ratio for Bs ^^/^^ (three dashed lines from left to right: 7 x 10^*^, 
2 X 10^^, 1 X 10^^) for tan/? = 50 in the mo-mi^2 plane. Other mSUGRA parameters are fixed to 
be ^0 = and fi> 0. The three short solid lines indicate the cr^o_p values (from left: 0.05 xlO^^ 
pb, 0.004 xlO^^ pb, 0.002 xlO"^ pb). The vertical dotted lines label Higgs masses. 
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FIG. 5. Branching ratio for ^^^^ (three dashed lines from left to right: 1.9 x 10^^, 

1 X 10^*^, 0.7 X 10^^) at tan /3 = 40 in the mo-m-x/2 plane. Other mSUGRA parameters are fixed to 
be ^0 = and jj, > 0. The three short solid lines indicate the cr^o_p values (from left: 0.03 xlO~^ 
pb, 0.002 xlO"^ pb, 0.001 xlO"*^ pb). The vertical dotted lines label Higgs masses. 
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FIG. 6. Branching ratio for Bg at tan/3 = 40 in the mQ-mi/2 plane for Aq = —2mi/2 

and > 0. The two short solid lines indicate the (7^o_p values (from left: 0.005x10"^ pb, 0.001 
xl0~^ pb). The vertical dotted lines label Higgs masses. 



12 



00 

o 
X 

T 

pq 




250 500 750 



1000 



FIG. 7. Branching ratio for Bg 
800 GeVin a ^ SUSY scenario (Ai22 



/i^/U as a function of mi/2 (in GeV) for mo = 300, 500, and 



A' 



i32 



0.02). Other mSUGRA parameters are fixed to be 



tan/3 = 10, ^0 = ^-nd /x > 0. Dashed hnes are to indicate the models that are excluded via the 
6 ^ S7 constraints. 
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